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Relativistic heavy ion collisions are studied assuming that particles can be described by a hadron 
gas in thermal and chemical equilibrium. The exact conservation of baryon number, strangeness 
and charge are explicitly taken into account. For heavy ions the effect arising from the neutron 
surplus becomes important and leads to a substantial increase in e.g. the n" /n^ ratio. A method 
is developed which is very well suited for the study of small systems. 



I. INTRODUCTION 



It has been shown recently in a very convincing manner that hadronic matter produced in coUisions of small 
T-H ' systems {e^e~ , p — p — p) is very close to chemical equilibrium. With only a small number of parameters, namely 
^ \ the temperature T, the volume V and, a factor 75 which measures the deviation of strange particles from chemical 
• equilibrium, it is possible to fit close to 30 different particle abundances . It is of great interest to establish whether 
or not this also holds for systems of medium size like S — S and for relativistic heavy ion collisions like Au — Au at 
the BNL-AGS or Pb — Pb at CERN This is of relevance because it is widely believed that in these collisions 
' one could create a quark-gluon plasma. Of immediate interest is the observed increase in the abundance of strange 
f — , particles in heavy ion collisions. In this paper we would like to present a method for taking into account the exact 
' conservation of baryon number B, strangeness S and charge Q. The method is ideally suited to analyze the extension 
_ from small systems like p — p towards bigger systems up to 5 — S. For larger systems the method becomes impractical 
-j-J • and should be replaced by a more straightforward approach based on the grand canonical ensemble. Our numerical 
,— H I results can be compared with data from the Brookhaven National Laboratory, e.g. the E866 collaboration [Q which 
^ , measures the dependence of hadronic ratios on the number of projectile participants thereby allowing one to study 
the transition from small to large systems in a systematic way. These results give insight into the behavior of the 
produced hadronic system as a function its size. The treatment presented in this paper differs from the standard 
one based on the grand canonical ensemble (see for example in that we consider the quantum number content 
exactly. This means that wc do not introduce chemical potentials for the baryon number or for strangeness. Chemical 
potentials are usually introduced to enforce the right quantum numbers of the system in an average sense. This is 
[ a correct treatment for large systems, however, for small systems the production of extra proton - anti-proton pairs 
for example will clearly be more suppressed than in large systems. These extra corrections were first pointed out by 
Hagedorn ||^ and subsequently a complete treatment has been formulated |p|~pj[| . We emphasize that these corrections 
do not contain information about the dynamics. They simply follow from baryon number conservation and must be 
taken into account before considering more involved models. It is also worth emphasizing that they do not introduce 
any new parameters. Our treatment differs from the one presented in in that it is much more analytic. It extends 
the analysis started in jl^ since it also includes the exact treatment of charge conservation and allows for the correct 
treatment of systems which are not isospin-symmetric. 



II. PARTITION FUNCTION 

The exact treatment of quantum numbers in statistical mechanics has been well established for some time |^ . It is 
obtained by projecting the partition function onto the desired values of B, Q and S 

Zb,q,s = ^ J^^ d(p e-^^* i- j^J da e"^^"^ fj di,e-^^^Z{T, Xb, XqAs) (2.1) 
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where Z is the (grand canonical) partition function and the usual fugacity factors As and As have been replaced by : 

\B=e}'^ XQ = e'" Xs^e'^. (2.2) 

As the contributions always come pairwise for particle and anti-particle, the fugacity factors will give rise to the cosine 
of the angle. In the extended treatment it is useful to group all particles appearing in the Particle Data Booklet 
into fourteen categories depending on their quantum numbers (we leave out charm and bottom). Zk is the sum (given 
below) of all mesons having strangeness ±1 {K, K, K*, . . .) and zero charge. Similarly Zn is the sum of all baryons 
and anti-baryons having zero strangeness, Zy is the sum of all hyperons and anti-hyperons, while Zq is the sum of all 
non-strange mesons. 
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Table 1 : Classification of the families of particles. 



For example, is the single particle partition function of all charged strange mesons {K^ , K , K*^, K* ...) i.e. 

je|s|=i,|B|=o,|Q|=i ^ ' 

where gj is the degeneracy factor. 

With the help of the definitions presented in table 1, the partition function can be rewritten in the following form, 

^ plTT />27r />27r 

Zb s q {T, V) — „ / d(f> exp (— i50) / dtp exp {—iBip) / da exp (—iQa) 
(27r)'^ Jo Jo Jq 

■ exp [2Zk cos (p + 2Zn cos tp + 2Z^c cos a + 2Z\ cos(?/; — <p) 

+2Zk- cos {(p + a) + 2Za- cos (i/j - a) + 2Zp cos {ip + a) + 2Za++ cos {ii + 2a) 

-|-2Zs+ cos {~(p + V + a) + cos {-(p + ip - a) 

+2Ze- cos {-2(p + ip-a) + 2Zho cos (-20 -I- ip) 

2Zacos(-3(/)-|-^-a)] ■exp(Zo). (2.4) 

To continue we would like to make use of the integral representation of the modified Bessel function 

/„(z) = i/ e^™'''' cos n6i.d6' (2.5) 
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However, this cannot be used in a straightforward manner since the dependence on angles does not factorize and 
almost all terms involve more than one angle. To circumvent this difficulty we introduce new angles whenever more 
than one appears. For example, for the term involving Z^ we introduce an intermediate angle A in the following way 
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(2.6) 
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With this method complete factorization can be achieved at the price of introducing a new summation at each step. 
Since there are ten cosines in Eq. (2.4) involving more than one angle, we end up having ten summations. The 
following result is obtained [| 
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(2.7) 



where 



vs ^ S — rii + n2 — 715 — Tie — '^nr — 2ng ~ Srig 

I'B = -B + ni + ^3 + n4 + 715 + ng + n7 + ng + ng + ?^io 

i^Q = Q + ^2 + «3 - + ^5 - '^e - '^s - '^g + 2nio. 



(2.8) 
(2.9) 
(2.10) 



Equation (2.7) is our main result and forms the basis of our analysis below. The sums in the above equation are 
dominated by terms corresponding to the index of the Bessel functions being zero (for the size of the system being 
sufhciently small) . This means however that for large values of the baryon number _B some of the indices could become 
quite large. In such cases evaluation of the sums becomes very time-consuming numerically and it is better to resort 
to the grand canonical ensemble. 

The canonical partition function for a gas with three conserved quantum numbers has thus been derived and will be 
used in the following to derive expressions for particle numbers. The differentiation of the equation (2.7) for particle 
abundances decreases/increases some of the indices vg^VB and vq by one and leads to useful R- factors which, for 
different particle species deviate from Zb.s,q^^ ^) in the way indicated below 
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If a particle, z, has strangeness 1, baryon number and charge 0, it's density will be given by 
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(2.13) 



^Full details can be found in [ p^ . 
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All other particle densities are obtained by using the appropriate R factor. The factor in square brackets in equations 
(2.12) and (2.13) replaces the fugacity in the usual grand canonical ensemble treatment Having thus determined 

all particle densities, we consider the behavior at freeze-out time. In this case all the resonances in the gas are allowed 
to decay into lighter stable particles. This means that each particle density is multiplied with its appropriate branching 
ratio (indicated by Br below). The abundances of particles in the final state are thus determined by : 

uh = n,Br{i H) (2.14) 

i 

where each sum runs over all particles contained in the hadronic gas and H refers to a hadron (7r+. , . . .). 



III. NUMERICAL RESULTS 



In order to stay close to the grand canonical ensemble treatment we keep the temperature and the density B/V 
fixed. The latter corresponds to keeping the baryon chemical potential fixed. We therefore expect all ratios to 
become constant as the size of the system is increased. Fig. 1 illustrates the approach to the thermodynamic (grand 
canonical) limit for the /'k'^ ratio. As the system becomes larger the ratio approaches a constant. The variation 
of the TT^/Tr"*" ratio as one increases the neutron surplus or, equivalently, the ratio B/2Q, is shown in Fig. 2. As 
expected one observes a clear increase of this ratio with B/2Q. In Fig. 3 we show the p/ir^ ratio as a function of the 
baryon number B. The /K~ ratio is shown in Fig. 4 for a fixed value of the temperature, T = 110 MeV. One 
observes a very smooth decrease and this ratio very quickly approaches its asymptotic value. 

Fig. 5 shows the p/tt"*" ratio as a function of the size of the system. This figure is interesting because preliminary 
results from the E866 collaboration Q| show a decrease of this ratio with increasing system size, i.e. the opposite 
behavior from the one observed in Fig. 5. In our opinion, this proves in a convincing manner that p are definitely not 
in thermal equilibrium. It would be difficult to show this if the analysis is based on only one system. 



IV. SUMMARY 



We have assumed that, as a first approximation, the hadronic final state produced in a relativistic heavy ion colli- 
sion can be described by a hadronic gas which is in chemical equilibrium. The standard description of such a system 
makes use of a statistical description based on the grand canonical ensemble using chemical potentials for the con- 
served quantum numbers. However, if the system is very small then there will always be corrections due to the size 
of the system. We have presented a method for taking into account corrections arising from the exact conservation of 
baryon number B, strangeness S and, charge Q. This method is well suited for small systems and numerical results 
have been presented for the most abundantly observed hadrons. For large systems the numerical evaluation becomes 
very slow because a large number of terms have to be kept in the sums. The numerical results presented can be 
compared with data from Brookhaven, e.g. the E866 collaboration M. We plan to do this in the near future. 
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FIG. 1. The /Tr~^ ratio as a function of the baryon number B for fixed values of the temperature T, the baryon 

density B/V and B/2Q. 

FIG. 2. The 7r~/7r+ ratio as a function of the B/2Q for a fixed value of the temperature T, the baryon density B/V, 

and the baryon number {B = 10). 

FIG. 3. The p/tt~^ ratio as a function of the baryon number B for fixed values of the temperature T, the baryon 

density B/V and B/2Q. 

FIG. 4. The / K~ ratio as a function of the the baryon number B for fixed values of the temperature T, the 

baryon density B/V and B/2Q. 

FIG. 5. The p/7r+ ratio as a function of the baryon number B for fixed values of the temperature T, the baryon 

density B/V and B/2Q. 



5 



0„ 24 



0„ 22 - 



□ „ 20 - 



□ „ 1 



□ „ 16 



0„ 1 



0„ 12 



„ 1 



T = 110 tle^ 
B/V = 0„ 05/fm 

B/2Q =1.5 



3 



2. 



4„ 



6„ 



8„ 



10„ 



B 




B/2Q 




B 




B 



0. 0002 



p /n 



^ ^ WO HeV 
B/V = 0„ 05/fm- 
B/2Q = 1 „ 
B/2Q = 1 „ 5 



0. 0001 



2.0 4.0 6.0^ 8.0 10. 



